Hole spin relaxation in [001] strained asymmetric Si/SiGe and Ge/SiGe quantum wells 
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Hole spin relaxation in [001] strained asymmetric Si/Sio.rGeo.a (Ge/Sio.aGeo.?) quantum wells 
is investigated in the situation with only the lowest hole subband being relevant. The effective 
Hamiltonian of the lowest hole subband is obtained by the subband Lowdin perturbation method in 
the framework of the six-band Luttinger k • p model, with sufficient basis functions included. The 
lowest hole subband in Si/SiGe quantum wells is light-hole like with the Rashba spin-orbit coupling 
term depending on momentum both linearly and cubically, while that in Ge/SiGe quantum wells is a 
heavy hole state with the Rashba spin-orbit coupling term depending on momentum only cubically. 
The hole spin relaxation is investigated by means of the fully microscopic kinetic spin Bloch equation 
approach, with all the relevant scatterings considered. It is found that the hole-phonon scattering is 
very weak, which makes the hole-hole Coulomb scattering become very important. The hole system 
in Si/SiGe quantum wells is generally in the strong scattering limit, while that in Ge/SiGe quantum 
wells can be in either the strong or the weak scattering limit. The Coulomb scattering leads to 
a peak in both the temperature and hole density dependences of spin relaxation time in Si/SiGe 
quantum wells, located around the crossover between the degenerate and nondegenerate regimes. 
Nevertheless, the Coulomb scattering leads to not only a peak but also a valley in the temperature 
dependence of spin relaxation time in Ge/SiGe quantum wells. The valley is actually due to the 
crossover from the weak to strong scattering limit. The hole-impurity scattering influences the 
spin relaxation effectively. In the strong (weak) scattering limit, the spin relaxation time increases 
(decreases) with increasing impurity density. The spin relaxation time is found to be on the order of 
1~100 ps (0.1~10 ps) in Si/SiGe (Ge/SiGe) quantum wells, for the temperatures, carrier/impurity 
densities and gate voltages of our consideration. 

PACS numbers: 72.25.Rb, 71.10.-w, 71.70.Ej 



I. INTRODUCTION 



In recent years, great efforts have been devoted to 
the design/realization of spintronic devices, which em- 
ploy the spin degree of freedom in traditional electron- 
ics for the sake of higher power efficiency, higher speed, 
and also greater functionality.—'^'^ Among different kinds 
of hosts for such devices, Si appears to be a particu- 
larly promising one and attracts much attention, partly 
due to the high possibility of eliminating hyperfine cou- 
plings by isotopic purification and the well developed 
microfabrication technology.^ In fact, the electron spin 
relaxation, which is necessary to be understood for the 
device design, has been widely investigated in Si ma- 
terials during the last decade. The study on relax- 
ation of electron spin qubit in Si quantum dot suggests 
that the relaxation rate can be strongly decreased by 
adding strain;^ The electron spin relaxation in asymmet- 
ric n-type Si/SiGe quantum wells (QWs) has been in- 
vestigated both theoretically^ and experimentally.--'^ It 
is shown that the electron spin relaxation time can be 
quite long (on the order of 10~^ ~ 10~^ s)r'2^ due to 
the weak Rashba spin-orbit coupling^ (typically about 
three orders of magnitude smaller than that in QW struc- 
tures based on III-V semiconductors -) . The electron spin 
transport/diffusion in bulk Si with a magnetic field per- 
pendicular to both the directions of spin polarization and 
spin transport /diffusion has also been studied recently. 



It is revealed that even in the absence of the traditional 
D'yakonov-Perel' (DP) relaxation mechanism,— there is 
an obvious spin relaxation along spin transport /diffusion, 
as predicted several years ago from a general QW model 
without any DP relaxation mechanism but with a mag- 
netic field in the Voigt configuration^^ That is also the 
case in the symmetric Si/SiGe QWsj^^ 

Although a broad interest has been taken in the 
electron spin relaxation in Si, to our knowledge, the 
hole spin relaxation has been rarely investigated so 
far. Glavin and Kim have calculated the spin relax- 
ation of two-dimensional holes in strained asymmet- 
ric Si/SiGe (Ge/SiGe) QWs four years agOfi^ and ob- 
tained a spin relaxation time of several tens of picosec- 
onds (several sub-picoseconds) in Si/SiGe (Ge/SiGe) 
QWs with large gate voltage (which induces an elec- 
tric field at 50 ~ 500 kV/cm) at room temperature. 
However, the results were obtained by means of the 
single-particle approximation,-'^^ therefore the effect of 
the carrier-carrier Coulomb scattering on spin relax- 
ation, which has been revealed to be important in spin 
relaxation )^^i^'^i^^i^^i^°i^^i^^ was not included. Besides, 
the nondegenerate perturbation method with only the 
lowest unperturbed subband of each hole state consid- 
ered as basis function is utilized to calculate the sub- 
band energy spectrum and envelope functions in Ref. [3 . 
However, as shown later in this paper, only consider- 
ing the lowest unperturbed subband is inadequate in 
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converging the calculation, but when more unperturbed 
subbands are included as basis functions, the nondegen- 
erate perturbation method fails. This work is to per- 
form a detailed investigation on hole spin relaxation in 
asymmetric Si/SiGe and Ge/SiGc QWs by means of the 
fully microscopic kinetic spin Bloch equation (KSBE) 
approach, with all the relevant scatterings included. 
Meanwhile, we apply the exact diagonalization method 
to obtain the energy spectrum and envelope functions, 
with sufficient unperturbed subbands included. In the 
KSBE approach, the momentum-dependent spin preces- 
sions give rise to the inhomogencous broadening, with 
which any scattering (including the Coulomb scattering) 
leads to an irreversible spin relaxationp^^ This approach 
has been successfully applied to study spin dynamics in 
quantum wire,^^'^'^ QW^2'5^>i^i2^i^i21 and bulk^* semi- 
conductor structures. The current work reveals that the 
Coulomb scattering plays a much more important role 
in hole spin relaxation in Si/SiGe (Ge/SiGe) QWs. It 
leads to a peak in both the temperature and density 
dependences of spin relaxation time in Si/SiGe QWs, 
where holes are generally in the strong scattering limit. 
Nonetheless, it leads to not only a peak but also a valley 
in the temperature dependence of spin relaxation time 
in Ge/SiGe QWs, where with the change of temperature 
the holes in Ge/SiGe QWs can be in the either strong or 
weak scattering limit. Besides, the spin relaxation time 
can be effectively influenced by the hole-impurity scat- 
tering, which tends to weaken the effect of the Coulomb 
scattering mentioned above with the increase of impurity 
density. 

This paper is organized as follows. In Sec. II the effec- 
tive Hamiltonian of the lowest hole subband (we focus on 
the situations with only the lowest subband being rele- 
vant) in asymmetric Si/SiGe (Ge/SiGe) QWs is derived. 
In Sec. Ill the KSBEs are constructed and the hole spin 
relaxation in Si/SiGe (Ge/SiGe) QWs is investigated. Fi- 
nally, we conclude in Sec. IV. 



II. EFFECTIVE HAMILTONIAN 

We start our investigation from the p-type 
Si02/Si/Sio.7Geo.3 (SiOs/Ge/Sio.sGeo.r) QWs. The 
Si02/Si/Sio.7Geo.3 and SiOa/Ge/Sio.sGco.y QW struc- 
tures are illustrated in Fig. [TJ The Si (Ge) layer is 
[001] I |z grown with a wide width (> 10 nm). The Si02 
layer is assumed to be an infinite potential barrier. 
With the valence band discontinuity at the Si/Sio.rGeo.a 
(Ge/Sio.sGeo.r) interface 55 meV (« 200 meV)}i 
ignored due to the large gate voltage (inducing an 
electric field > 50 kV/cm) and wide well width^ii. the 
triangular potential approximation is adopte d^'*'^^i'^° 
and the well width then becomes irrelevant. 

Based on the theory of Luttinger-Kohn^ii^^ and Bir- 
Pikusr^ the valence-band structure of the strained 
QWs can be described by the 6x6 effective-mass 



TABLE I: Material parameters of Si and Ge. The mass den- 
sity d, deformation potentials Dac and Aop, optical phonon 
energy luoop and sound velocity Va are taken from Ref. Issl . 
The Luttinger parameters 71, 72, 73 are from Ref. [30I. 
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Hamiltonian fiii^i 
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Here = + H)]'' is the Luttinger 

Hamiltonian^iiHiH with H^^ corresponding to the 
part with k^^y = 0- is the contribution due to the 
biaxial strain.*^— I'^^i^'^i'^^ ^(^) is the confining potential 
and /g is the 6x6 unit matrix. The 2-components of the 
subband envelope functions (the x- and y-components 
are plane waves) obtained by solving the eigen-equation 
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Here Z, h and s represent the light hole (LH), heavy hole 
(HH), and spht-off (SO) hole states, respectively, and n 
is the subband number. The solution of the envelope 
functions is stated in Appendix A. 

The Lowdin partitioning'^- is performed upto second 
order in i/j^"^ on the basis constructed by ^Aan (A=l, 
2; a—h, I, s) to obtain the effective Hamiltonian of the 
lowest hole subband.^^ Due to the biaxial strain,i^i22i^i^ 
the lowest subband in Si/SiGe QWs is a LH-like state 
(LHO), which is an admixture of LH and SO hole states, 



3 



with 



e 



. Si02 


Si 


SiQ7GeQ3 




1 i 

i i j> 

i 1 
/ "\ 

(■ \ 


(a) 








t 


Ge 

•' > ^"-"^ 






! 1 

! 1 
! 1 

1 

■ \ 











0.4 




3 




0.2 




0.1 


3 
.c 




w 









■s 




a 

1 


0.3 


1 






0.2 








0.1 










-10 



10 15 20 25 30 35 
z (nm) 



FIG. 1: (Color online) Schematics of the Si02/Si/Sio.7Geo.3 
QW structure (a) and Si02/Ge/Sio.3Geo.7 QW structure (b). 
Two vertical dashed lines in each figure represent the two in- 
terfaces. The solid curves represent the confining potential 
V{z) with electric field E = 300 kV/cm. The valence band 
discontinuities at the Si/Sio.rGeo.a and Ge/Sio.sGeo.r inter- 
faces are neglected in the triangular potential approximation. 
The chain curves with their scale on the right hand side of the 
frame are |*I'iio(2)P (a) and I^PihoC^)!^ (b), respectively, rep- 
resenting the lowest LH and HH distributions in Si/Sio.rGeo.a 
and Ge/Sio.3Geo.7 QWs along the z-axis. 
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while that in Ge/SiGe QWs is a pure HH state (HHO). 
The effective Hamiltonian of the lowest hole subband in 
Si/SiGe (Ge/SiGe) QWs can be written asii 
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where k is the in-plane momentum, m*^'-' [m is the 
in-plane effective mass of the lowest light (heavy) hole 
subband in Si/SiGe (Ge/SiGe) QWs, cr are the PauH 
matrices, and fl'-^^ [il^'*)] is the Rashba term of the LHO 
(HHO) subband in Si/SiGe (Ge/SiGe) QWs. fl^') has 
both the linear and cubic dependences on momentum, 
whereas fl^'^^ has only the cubic dependence. For the 
LHO subband in Si/SiGe QWs, 
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Here A, B and C are the valence band parameters, which 
relate to the Luttinger parameters (Table |T| 71 , 72 and 
73 through A = 71, B = 272 and V3B2 4, (j2 = 2\/373. 



{a—h, I, s) are the subband energy levels. X^^^'' 
kI"'^'' are defined as A^"'?^ = dzxi^^\z)xl^,\z) and 
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It is noted that in Ref . UA 



the coefficients H, Q and A miss the pre-factor — i and 
meanwhile H misses the summation over I (i.e., the con- 
tribution due to the higher LH subbands, which is in fact 
negligibly small). Unlike the work by Glavin and Kim— 
where the nondegenerate perturbation method with only 
the lowest unperturbed subband of each hole state be- 
ing accounted [refer to Eqs. (3-9) in Ref. T^] is employed 
in obtaining the subband energy spectrum E'i"'' and the 
envelope functions ^'aqm, we apply the exact diagonaliza- 
tion method with sufficient unperturbed subbands [2 (20) 
for each hole state in Si/SiGe (Ge/ SiGe) QWs] for the 
sake of convergence. In fact, the calculation with only the 
lowest unperturbed subband of each hole state is inade- 
quate for the convergence. Moreover, the calculation with 
more unperturbed subbands included may cause diver- 
gence as long as the nondegenerate perturbation method 
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is utilized. When performing the subband Lowdin par- 
tition method^ to obtain the effective Hamiltonian of 
the lowest subband, we choose sufficient envelope func- 
tions '^\an as basis functions. For Si/SiGe (Ge/SiGe) 
QWs, the number of envelope functions containing "^xin 
and is 4 (40) in total, and the number of envelope 

functions 'i'xhn is 2 (20). The reason that more basis 
functions are needed for Ge/SiGe QWs comes from the 
fact that the couplings between the hole subbands are 
stronger than those in the Si/SiGe QWs (the Luttinger 
parameters in Ge are much larger than those in Si but 
under the strain the energy differences between hole sub- 
bands are comparable in Si/SiGe and Ge/SiGe QWs). 
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FIG. 2: (Color online) (a) The subband energy spectrum at 
r point against electric field in Si/SiGe QWs. Four subbands 
are shown: the first (second) light-hole subband LHO (LHl), 
the first heavy-hole subband HHO and the first split-off hole 
subband SOO. (b) Spin-orbit coupling coefficients ^E,, ^11 
and |e for LHO subband in Si/SiGe QWs and |A for HHO 
subband in Ge/SiGe QWs against the electric field E. The 
scale of |H is on the right hand side of the frame. 

According to Eqs. (|3|) and (|12p. the in- plane effective 
mass of LHO subband in Si/SiGe QWs m''-' is calculated 
to be about 0.27too in the whole electric field range under 
consideration, and that of the HHO subband in Ge/SiGe 
QWs m^''^ is 0.057mo. We plot the energy levels of four 
subbands in Si/SiGe QWs (the first and second LH sub- 
bands LHO and LHl, the first HH subband HHO, and 
the first SO subband SOO) at F point in Fig. ^a.) and 
the spin-orbit coupling coefficients of the LHO subband 



in Si/SiGe QWs (S, H and 9) and the HHO subband in 
Ge/SiGe QWs (A) in Fig.^b). As shown in Fig.^a), a 
crossing between the LHl and HHO subbands appears at 
about E = 360 kV/cm and an anticrossing between the 
LHl and SOO subbands appears around E ~ 500 kV/cm, 
respectively. It is noted that notwithstanding the fact 
that our calculation goes to the infinitesimal electric field 
regime in Fig. [21 only the results in the large electric field 
regime (i.e., E' > 50 kV/cm for QWs with well width 
~ 10 nm) are valid (the spin relaxation investigated later 
is also in the large electric field regime), as our model 
fails in the small electric field regime due to the disre- 
gard of the discontinuity in the Si/SiGe (Ge/SiGe) inter- 
face. In this work the temperature dependence of band 
parameters^i^ is not taken into account, due to both the 
weak temperature dependence of band parameters and 
the negligible contribution from the conduction band.— 
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FIG. 3: (Color online) (a) Spin splitting for holes with in- 
plane kinetic energy being fcsT (T = 300 K) along the [100] 
(the solid and dotted curves) and [110] (the chain and dashed 
curves) directions in Si/SiGe QWs. Solid and chain curves: 
results of our calculation; dotted and dashed curves: results 
from Ref. (b) Spin-orbit coupling coefficient |n calcu- 
lated with the nondegenerate perturbation and exact diago- 
nalization methods respectively. The solid (chain) curve and 
the dots (squares) are results from the perturbation and ex- 
act diagonalization methods, respectively with the lowest one 
(two) unperturbed subband (subbands) of each hole state con- 
sidered. 

We also calculate the spin splitting for holes in Si/SiGe 
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QWs with in-plane kinetic energy being fc^T (T = 
300 K) along the [100] (solid curve) and [1101 (chain 
curve) direction in Fig. El^a), as done in Ref. Ilj- The 
dotted and dashed curves are taken from Ref. ^M', corre- 
sponding to the spin splittings along the [100] and [110] 
directions, respectively. It is shown that our results differ 
from those in Ref. IT^4^ 

We examine our results by further carrying out a cal- 
culation of spin-orbit coupling coefScients/spin splitting 
with the envelope functions ^xan obtained by the nonde- 
generate perturbation method as that in Ref. JLi (also re- 
fer to Appendix A of this paper), but with the spin-orbit 
couphng coefficients obtained in this work [i.e., Eqs. ^ 
[TT|l and P^ ]. As a comparison, we plot the electric field 
dependence of 11 [Eq. (fTO|l ] in Fig. [31(b), where the solid 
(chain) curve and the dots (squares) are results from the 
perturbation and exact diagonalization methods, respec- 
tively with the lowest one (two) unperturbed subband 
(subbands) of each hole state considered. We find that 
when only the lowest unperturbed subband of each hole 
state is considered, the perturbation calculation and our 
exact diagonalization calculation yield almost the iden- 
tical results [compare the solid curve and the dots in 
Fig. [IKb)]. We also find that, as said above, when more 
unperturbed subbands of each hole state are accounted, 
the nondegenerate perturbation method may cause diver- 
gence in the spin-orbit coupling coefficients/spin splitting 
(a divergence near E — 500 kV/cm in the chain curve is 
observed). The divergence is caused by the degeneracy 
of the LH subband and the SO subband (refer to Ap- 
pendix A for details), and disappears in the exact diago- 
nalization calculation (see squares in the figure). Besides, 
the large discrepancy between the results with different 
number of unperturbed subbands included indicates that 
only considering the lowest subband of each hole state is 
inadequate for the convergence of calculation. As a re- 
sult, the exact diagonalization calculation with sufficient 
unperturbed subbands of each hole state included is nec- 
essary. 



III. HOLE SPIN RELAXATION 

We perform the fully microscopic KSBE approaches to 
study the hole spin relaxation. The KSBEs constructed 
by the nonequilibrium Green function method read^^ 



d , . d , , 



coh 



(17) 



in which pk represent the density matrices of holes. 
^Pk(i)|j,Q[j are the coherent terms describing the coher- 
ent spin precessions due to the effective magnetic fields 
from the Rashba term and the Hartree-Fock Coulomb in- 
teraction. ^Pk(^)|scat ^tand for the scattering terms, in- 
cluding the hole-deformation optical/acoustic phonon;^ 
hole-impurity and hole-hole Coulomb scatterings. Ex- 
pressions of the coherent and scattering terms are given 



in detail in Ref. '27'. What need to be specified are the ma- 
trix elements of the hole-phonon interaction in the scat- 
tering terms. The matrix elements of hole-deformation 
acoustic phonon scattering and hole-deformation opti- 
cal phonon scattering are IMqcqI^ = ^2dv^ [-^(^g^)!^ and 



IM, 



op.Ql 



2duJo 



|/(igz)|^ respectively. Here Q = (q, g^) 
is the phonon momentum. The values of mass density d, 
deformation potentials Dac and Aop, the optical phonon 
energy hujop and the sound velocity Vg in Si and Ge are 
listed in TablelH |/(«g^)|^ is the form factor with I{iqz) — 

/-!l*U(^)e''^^*i/o(^)rf^ - {x'i'\^)W'''\xi'\z)) + 
{xV{z)W'^'"\xViz)) for Si/SiCe QWs and I{iq,) ^ 
!^^^U{z)e^'^^'^im{z)dz = {x'^\z)W'^'^\x^^\z)) for 
Ge/SiCe QWs. By numerically solving the KSBEs, one 
can obtain the time evolution of density matrices and 
then the spin relaxation time. In the calculation, the 
initial spin polarization of holes is set to be 5 %. 



A. Hole spin relaxation in Si/SiGe QWs 

We first study the spin relaxation of the lowest hole 
subband in Si/Sio.7Geo.3 QWs. E = 300 kV/cm un- 
less otherwise specified. The LHO holes have a distri- 
bution along the z-direction as shown in Fig. [UJa) by 
the chain curve. The spin-orbit coupling coefficients 
are |S = 1.33 meV nm, |n = -0.83 meV nm^ and 
fe = -0.39 meV nm^, respectively [Fig. [^Jb)]. More- 
over, due to the quite small material parameter B (or 
72) in Si, the linear part of the Rashba term in Si/SiGe 
QWs is relatively more important. The main results are 
plotted in Figs.HHTl showing the spin relaxation with dif- 
ferent temperatures, carrier/impurity densities and scat- 
terings. 
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FIG. 4: (Color online) Spin relaxation time r against temper- 
ature T with different hole densities. The impurity density 
Ni = Q and the electric field E = 300 kV/cm. 

Figure m shows the temperature dependence of spin re- 
laxation time with different hole densities. The impurity 
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density is set to be zero. A peak, appearing at a tem- 
perature around the Fermi temperature Tj* = E^/ks 
[E'j is the hole Fermi energy. Tj* w 35 K with den- 
sity Nh = 4 X 10^"'^ cm~^), is observed except when the 
hole density is too low. This kind of peak has been pre- 
dicted by Zhou et al. in high-mobility n-doped GaAs 
QWgis and later observed by Ruan et al. experimentally 
at about Tj/2 in the temperature dependence of electron 
spin relaxationji^ Similar peaks have also been predicted 
very recently in the temperature dependence of electron 
spin relaxation at a temperature in the range of (Tj/4, 
Tj/2) in intrinsic bulk GaAs2^ and at a temperature 
around the hole Fermi temperature in impurity-free 
p-type GaAs QWs where the hole density is much higher 
than the electron densityj^ In fact, this feature ap- 
pears in the electron spin relaxation of strong scattering 
system with the DP relaxation mechanism being domi- 
nant when the Coulomb scattering (either the intraband 
electron-electron or the interband electron-hole Coulomb 
scattering) is the main scattering. When the intraband 
electron-electron Coulomb scattering dominates ) 
the peak appears around the crossover from the degen- 
erate to nondegenerate regime of electrons. When the 
interband Coulomb scattering dominates, the peak ap- 
pears around the crossover from the degenerate to non- 
degenerate regime of holes in p-type systemsi^ 

It is known that in the strong scattering system, 
strengthening scattering can suppress the inhomogeneous 
broadening and tends to prolong the spin relaxation time 
within the DP relaxation mechanism )^^i^^i^^i^^'^'^i^^ and 
that the Coulomb scattering rate has a dependence in 
the degenerate regime but a (T"'^/^) dependence in 
the nondegenerate regime in the two (three)-dimensional 
carrier systemsJ^i^ Thus with the increase of T, the 
dominant Coulomb scattering tends to cause first an in- 
crease and then a decrease in the electron spin relaxation 
time. Meanwhile, the increase of inhomogeneous broad- 
ening with T tends to cause a monotonous decrease in 
the spin relaxation time and thus a shift of the peak in 
the T-T curve towards the lower temperature. The mag- 
nitude of the latter effect depends on the form of the 
momentum dependence of the DP term. When the DP 
term mainly depends on the momentum linearly (cubi- 
cally), the latter effect is moderate (strong). The above 
scenario also holds in the hole spin relaxation, such as the 
case considered here. The hole system in Si/SiGe QWs 
is in the strong scattering limit where the Coulomb scat- 
tering dominates (as discussed later) and the linear part 
of the Rashba term is more important, thus the peak in 
the T-T curve is obvious near T^. However, when the 
hole density is low enough and thus holes are in the non- 
degenerate regime throughout the temperature regime 
under consideration, the spin relaxation time decreases 
monotonously with temperature, as shown by the solid 
curve in Fig. U) 

A similar phenomenon is expected to happen in the 
density dependence of spin relaxation time, as the 
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FIG. 5: (Color online) Spin relaxation time r against hole 
density Nh with different temperatures. The impurity density 
Ni = and the electric field E = 300 kV/cm. 



Coulomb scattering rate has an Nf^-^ i^h^^^) depen- 
dence in the degenerate regime while an Nh dependence 
in the nondegenerate regime in two (three)-dimensional 
systemsiiii^ This is exactly the case ;^^'^^ as shown by 
Fig. [5l Nevertheless, when the temperate is high enough 
and thus holes are in the nondegenerate regime through- 
out the density range under consideration, the spin relax- 
ation time increases monotonously with density, as shown 
by the dotted curve in Fig. [5l 

In Fig. [6] the spin relaxation time against temperature 
with different impurity densities is plotted. It shows that 
adding impurities reduces spin relaxation rate. That is 
because the inhomogeneous broadening is suppressed by 
introducing hole-impurity scattering in the strong scat- 
tering limit With the increase of impurity 
density Ni, the hole- impurity scattering, which is insensi- 
tive to T in low temperature regime, becomes important. 
Thus the peak in r-T curve due to the Coulomb scat- 
tering becomes less pronounced or even disappearsJ^i^^ 
That is the reason why the peak is easier to be observed 
experimentally in high-mobility samplesi ^^'^^ It is also 
noted that the impurity scattering has marginal effect 
on spin relaxation near room temperature, which is un- 
derstood by recalling that the impurity scattering rate 
decreases with carrier energy4^ 

To understand the relative importance of different 
scatterings in spin relaxation, we calculate the spin relax- 
ation time with different scatterings included and show 
its density dependence in Fig. [71 T is taken to be 300 K. 
In Fig. [TKa) the impurity density Ni = 0. The solid 
curve corresponds to the case with all the scatterings 
(the hole-hole Coulomb, hole-optical phonon and hole- 
acoustic phonon scatterings) included. The dashed, dot- 
ted and chain curves correspond to the cases without the 
hole-optical phonon, hole-acoustic phonon and hole-hole 
Coulomb scattering, respectively. By comparing these 
four curves, one finds that: (i) Even with T = 300 K, 
the hole-hole Coulomb scattering plays a much more im- 
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FIG. 6: (Color online) Spin relaxation time r against temper- 
ature T with different impurity densities. The hole density 
Nh = 4:X 10" cm-2 and the electric field E = 300 kV/cm. 



portant role than the hole-phonon scattering (in addi- 
tion, similar calculations show that when T < 200 K, 
the hole-phonon scattering can be completely ignored); 
(ii) The acoustic phonon scattering plays a relatively 
more efficient role than the optical phonon scattering 
[that is because the optical phonon energy hujop is high 
(63 meV) while the hole Fermi energy is low due to the 
large in-plane effective mass (0.27mo)]. In fact, the hole 
system under consideration is in the strong scattering 
limit generally, but falls into the weak scattering limit 
when the hole-hole Coulomb scattering is removed arti- 
ficially. (When T — 300 K here, the momentum relax- 
ation time Tp due to the relatively stronger hole-acoustic 
phonon scattering is about 2.1 ps, while the mean spin 
precession rateiSii^i^^ (fi^'^) is about 0.52^0.56 ps~^ with 
the change of hole density. Thus without the hole-hole 
Coulomb scattering, rp(ri('^) > 1, indicating the weak 
scattering limit.) This feature can be justified by com- 
paring two groups of curves in Fig.[7]Jb). One group with 
the hole-hole Coulomb scattering (the solid and dashed 
curves) indicates that adding a small amount of impuri- 
ties helps increasing r, corresponding to the strong scat- 
tering case, while the other group without the hole-hole 
Coulomb scattering (the chain and dotted curves) shows 
an inverse effect (i.e., a decrease in r) with adding quite 
a small amount of impurities, which typically happens in 
the weak scattering limit. 

Finally, we investigate the electric field dependence of 
spin relaxation. The hole spin relaxation time r against 
electric field E with different temperatures is plotted in 
Fig. [H For each temperature, we choose the appropriate 
range of electric field to ensure that the effect of sec- 
ond hole subband is irrelevant. It is shown that with 
the increase of electric field, the spin relaxation time is 
shortened due to the strengthened spin-orbit coupling. 
Besides, as a comparison to the features of spin relax- 
ation with E — 300 kV/cm, we also present the hole 
density and temperature dependences of spin relaxation 
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FIG. 7: (Color online) Spin relaxation time r against hole 
density Nh with different scatterings included. T — 300 K 
and E — 300 kV/cm. (a) Solid curve: with all scatter- 
ings [the hole-hole Coulomb (h-h), the hole-optical phonon 
(h-op) and the hole-acoustic phonon (h-ac) scatterings] in- 
cluded; Dashed curve: without the h-op scattering; Dotted 
curve: without the h-ac scattering; Chain curve: without the 
h-h scattering. All the four cases are calculated with impurity 
density Ni — 0. (b) Solid curve: with the h-h, h-op and h- 
ac scatterings; Dashed curve: same as the solid curve with 
the additional hole-impurity scattering {Ni = lO^'' cm~^) 
included; Chain curve: with the h-op and h-ac scatterings; 
Dotted curve: same as the chain curve with the additional 
hole-impurity scattering {Ni — 5 x 10* cm~^) included. 



with E ~ 600 kV/cm in Fig. [51 One finds that the peak 
in hole density/temperature dependence of the hole spin 
relaxation time (the dashed/sohd curve in Fig. [9]) due 
to the Coulomb scattering still exists. Moreover, the lo- 
cation of the peak remains almost the same despite the 
change of the gate voltage. This indicates that when the 
linear part of the Rashba term is important, the trend 
of variation of the hole spin relaxation time mainly as- 
sociates with that of the Coulomb scattering strength 
around the crossover between the degenerate and nonde- 
generate regimes, whereas the increase of the inhomoge- 
neous broadening with increasing temperature/hole den- 
sity is moderate (even with the increase of the spin-orbit 
coupling coefficients by the larger gate voltage). 
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FIG. 8; (Color online) Spin relaxation time r against electric 
field E with difi'erent temperatures. The hole density is Nh = 
10^^ cm~^ and the impurity density Ni = 0. 
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FIG. 9: (Color online) Spin relaxation with electric field E — 
600kV/cm. Dashed curve: spin relaxation time r against hole 
density A^^ with temperature T = 50 K; Solid curve: spin 
relaxation time r against temperature T with hole density 
Nh = 4 X 10^^ cm~^ (Note the scale is on the top of the 
frame). Ni=0. 



B. Spin relaxation in Ge/SiGe QWs 



relatively high, not only a peak in r-T curve is present, 
but also a valley before the peak is observed. 
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FIG. 10: (Color online) Spin relaxation time r against 
temperature T in Ge/SiGe QWs. The electric field E = 
300 kV/cm. The curve labeled by (a, 6) corresponds to the 
case with hole density Nh = ax 10^^ cm~^ and impurity den- 
sity Ni=bx 10" cm"^. The electric field is E = 300 kV/cm. 
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FIG. 11: (Color online) Spin relaxation time r against electric 
field E at different temperatures. The hole density is Nh = 
10^^ cm~^ and the impurity density Ni is zero. 



The hole spin relaxation in Ge/Sio.aGeo.y QWs is 
also investigated. The HHO holes have a distribution 
along the z-direction as shown in Fig. [Ijb) by the chain 
curve when the electric field is E — 300 kV/cm. In 
Fig. [TO] the spin relaxation time r against tempera- 
ture T with electric field E = 300 kV/cm [at which 
|A = —6.06 meV nm'^, as shown in Fig. [2Jb)] is plotted. 
It is shown that the hole spin relaxation time in Ge/SiGe 
QWs is much shorter than that in Si/SiGe QWs. That is 
because the inhomogeneous broadening in Ge/SiGe QWs 
is quite strong, as the spin-orbit interaction in Ge/SiGe 
QWs is relatively stronger (in consistence with the heav- 
ier Ge element) and the Rashba term depends on mo- 
mentum cubically. Apart from the fast spin relaxation, 
a new phenomenon emerges - when the hole density is 



Unlike the case in Si/SiGe QWs where the peak is close 
to the Fermi temperature (Fig. [4]), the peak here is lo- 
cated at a temperature about T^/2 (e.g., 170 K 
with density Nh = 4 x 10^^ cm~^). As discussed pre- 
viously about the spin relaxation in Si/SiGe QWs, the 
Coulomb scattering strength first increases and then de- 
creases with increasing temperature (accompanying the 
crossover from the degenerate to nondegenerate regime), 
tending to cause a peak in the r-T curve around the 
Fermi temperature in the strong scattering limit. How- 
ever, the enhancement of inhomogeneous broadening 
with T is strong here as the Rashba term depends on 
momentum cubically (unlike the case in Si/SiGe QWs 
where the linear part of the Rashba term is important). 
Thus with both effects accounted, the shift of the peak 
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towards a lower temperature is expected. By compar- 
ing the spin relaxation in high temperature regime of 
the three curves with A^^ = in Fig. [101 one finds that 
the spin relaxation time r decreases with hole density Nh 
monotonically, which is different from the case in Si/SiGe 
QWs (indicated by the dotted curve in Fig[5]). That is be- 
cause the Rashba term here depends on the momentum 
cubically and the inhomogeneous broadening increases 
with density strongly with an A^^ dependence. 

While the peak is associated with the crossover from 
the degenerate to nondegenerate regime, the valley in 
Fig. [TO] is actually related to the crossover from the weak 
to strong scattering limit. When the temperature is low 
enough and the density is high, the hole system is highly 
degenerate (due to the high Fermi energy) and thus the 
dominant Coulomb scattering becomes very weak. The 
hole system then falls into the weak scattering limit. 
Therefore with the increase of temperature from the very 
low temperature, the strengthening of scattering reduces 
the spin relaxation time first, until the crossover to the 
strong scattering system, and then increases the spin re- 
laxation time - this leads to a valley in the r-T curve 
in the degenerate regime. The dotted curve in Fig. [TO] 
stands for the case with the same hole density as the 
chain curve but with a small amount of impurities. It 
is observed again that introducing a weak impurity scat- 
tering in the weak (strong) scattering limit leads to a 
decrease (an increase) of the spin relaxation timeJ^^ 

The electric field dependence of spin relaxation is also 
investigated, with the hole spin relaxation time r against 
electric field E under different temperatures plotted in 
Fig. [TTJ It is shown that the spin relaxation time de- 
creases with electric field slowly in the large electric 
field regime, corresponding to the marginal electric field 
dependence of the spin-orbit coupling strength in the 
large electric field regime [refer to the dotted curve in 
Fig. [lb)]. 



IV. CONCLUSION 

In conclusion, we have investigated the hole spin 
relaxation in [001] strained asymmetric Si/Sio.yGeo.a 
(Ge/Sio.sGeo.r) QWs with large gate voltage in this work. 
We focus on the situations with only the lowest hole sub- 
band being relevant. The effective Hamiltonian of the 
lowest hole subband is obtained by the subband Lowdin 
perturbation method in the framework of the six-band 
Luttinger k • p model, with sufficient basis functions in- 
cluded for the convergence of calculation. Due to the 
biaxial strain, the lowest subband in Si/SiGe QWs is a 
light hole-like state, while that in Ge/SiGe QWs is a 
heavy hole state. 

We apply the fully microscopic KSBE approach to in- 
vestigate the hole spin relaxation in Si/SiGe (Ge/SiGe) 
QWs. By means of this approach, all the relevant scat- 
terings, such as the hole-phonon, hole-impurity and the 
hole-hole Coulomb scatterings can be taken into account 



explicitly. It is discovered that the hole-phonon scat- 
tering is very weak compared to the hole-hole Coulomb 
scattering, even at high temperatures. This makes the 
hole-hole Coulomb scattering to play a very important 
role in spin relaxation. It leads to a peak of spin relax- 
ation time in both the temperature and carrier density 
dependences in Si/SiGe QWs. The peak is associated 
with the crossover from the degenerate to nondegener- 
ate regime of hole system, and thus locates around the 
crossover point. However, the increase of inhomogeneous 
broadening with temperature/hole density tends to lead 
to a shift of the peak towards a lower temperature/hole 
density. The magnitude of the shift depends on the 
form of the momentum dependence of the Rashba term. 
For Si/SiGe (Ge/SiGe) QWs, the Rashba term mainly 
(only) depends on momentum linearly (cubically), and 
thus the shift of the peak is marginal (obvious). In ad- 
dition, in contrast to the GaAs QWs where the peak 
in the temperature dependence of the electron spin re- 
laxation can only be observed for high mobility sam- 
ples with low carrier densityji^ the peak predicted in 
Si/SiGe (Ge/SiGe) QWs can be observed even at high 
carrier density, thanks to the weak hole-phonon scatter- 
ing. The Coulomb scattering also leads to a valley at 
low temperature in the temperature dependence of hole 
spin relaxation time in Ge/SiGe QWs with high hole den- 
sity, which is related to the crossover from the weak to 
strong scattering limit. The hole spin relaxation time can 
be effectively influenced by the hole-impurity scattering, 
tending to weaken the effect of Coulomb scattering men- 
tioned above with the increase of impurity density. Apart 
from the abundant temperature and hole/impurity den- 
sity dependences, the spin relaxation time decreases with 
the gate voltage, accompanying the increase of spin-orbit 
coupling strength. 

The hole spin relaxation time, depending on the tem- 
perature, carrier/impurity density and gate voltage, is 
found to be on the order of 1~100 ps (0.1^10 ps) in 
Si/SiGe (Ge/SiGe) QWs within the scope of our inves- 
tigation. These time scales are much shorter than the 
electron spin relaxation time in Si/SiGe QWs (on the or- 
der of 10^^ ~ 10~^ s)i^iii^ A hole spin relaxation time 
on the order of 0.1~1 ps was theoretically reported in p- 
doped GaAs QWs (with temperature being 100^300 K, 
hole density being 0.5 ~ 4.5 x 10^^ cm^^, impurity den- 
sity being 0~1 times hole density and gate voltage in- 
duced electric field being about 100 kV/cm)^ and a hole 
spin relaxation time of 4 ps was experimentally observed 
in n-doped GaAs QWs (at 10 K)."*^ Thus, generally, the 
hole spin relaxation time in Si/SiGe (Ge/SiGe) QWs is 
longer than (comparable with) the hole spin relaxation 
time in GaAs QWs. It should be pointed out at last that 
the strain in Si/SiGe (Ge/SiGe) QWs plays an impor- 
tant role in spin relaxation, as it shifts the energy levels 
of the light hole states away from the heavy hole ones. 
If the strain is removed (e.g., in Si02/Si or Si02/Ge in- 
version layer), the coupling between the light hole and 
heavy holea states are strengthened and the hole spin 
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relaxation should be enhanced. 
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APPENDIX A: SOLUTION OF THE ENVELOPE 
FUNCTIONS *A«n 

The envelope functions ^'Acm(^) {^—^j 2; a—h, I, s) 
satisfy the eigen-equation 

Ho-^xaniz) = El^^^fxMz) (Al) 

with Ho = + i?e + V{z)Ie. ^'A/m(z) has only one 
component Xn \ satisfying [-^£7 + V{z)]xn\z) = 

Ei!^^Xji\z) with mi''^ = mQ/{A—B), which can be solved 
directly. The envelope functions ^\()n{z) {/3 =1, s) have 
two components, with Xri^"^ and Xn^"^^ satisfying 

(A2) 



(Is) 






Here 



^^^^ 2mo dz^ 2mo dz^ ^ 



E, 



V2 




V{z)h 



(A3) 



with the second term representing the contribution from 
the biaxial strain,^iHi33^ = -2&(2ci2/cii + 1)5, 
where cn and C12 are the elastic constants, b is the defor- 
mation potential constant and 5 is the relative lattice mis- 
match in the interfacci^i^ E^ is 95.8 meV (—115.7 meV) 
for Si/Sio.7Geo.3 (Ge/Sio.aGeo.r) QWs. A is the SO 
splitting. The solutions with {x!n^\z)\x\?^\z)) > 

{Xri'^\z)\Xn'^\z)) correspond to the LH-like states and 
thus (3 equals to I. Otherwise the solutions are deemed 
as the SO-like states with (3 = s. 

An unitary transformation which diagonalizes the 
strain term in Eq. (|A3|) is performed on Eq. (|A2[) , leading 



to Hq^^ as 



I 2mi dz^ 

I 

\ 2m* dz 

+ V{z)l2. 



E, 



2m* dz^ 
2m2 dz'^ 



(A4) 



U 



1 

/TvT 



W, E. I \ No E. 



1 

/7h 



(A5) 



is the unitary matrix, with 
A + ^9E^/4 + AE, + A^), E2 



El 



^9E^/4: + AE, + A^) and 7Vi,2 = 1 
mi, m2 and m* in Eq. (jA4[) are 



- A - 

2(1 - Ei^2/ES^ 



mi 



m2 



mo 



mo 



A + B - 



A + B I 



9EJ4 + A/2 
2 ' y/9E^/4 + AE, + A^ 

9EJA + A/2 



2 y^9E^/4 + AE, + A^ 



mo 



E,^9/A + A/E, + {A/E,)^ 
V2AB 



(A6) 

-1 

; 

(A7) 
(A8) 



It is noted that expressions (jAGIIASP about the effective 
masses are valid for both the Si/SiGe QWs (with E^ > 0) 
and Ge/SiGe QWs (with E,<0), while those in Ref.H 
[Eqs. (5-6)] equal to Eqs. (|A6IIA8p only when E^ > 0, i.e., 
they are valid only for Si/SiGe QWs. 



Hg''' can be separated into the diagonal Hg''j^' and off- 
diagonal Hqq parts. In the nondegenerate perturbation 

~ (Is) 

method, H^q' is treated as the perturbation term. The 

eigen-equation of iJoc' '^^^ solved directly, with the 
eigen-values Eni and En2 and the corresponding eigen- 
functions 



(Is) 




and 




(A9) 



determined by the Schrodinger equation [^^j^ + E^ + 

^{z)]xnd^) = En^Xniiz) (C = 1, 2). In this paper, 
the term "unperturbed subbands" mentioned in the dis- 
cussion of perturbation (exact diagonalization) method 

actually means the eigen-states of , i.e., the states 
with energy levels being -E„i,2 and wavefunctions being 
Eq. (1X91) . 

To the first order of Hqq\ the perturbed eigen-values 

are En"^ = E^i and E^^ = En2, and the corresponding 
perturbed eigen-functions are 



and 



Xnl{z) 



E„' -£?!^,xn'i{z) 

Xn2{z) 



(AlO) 



(All) 
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respectively. Here = ^zxn^-^ 

^' = 1, 2). Finally the eigen- values si^^ of H^^"^ 
are En^ and En \ and the corresponding cigen-functions 

(A12) 

and 

U I ) . (A13) 

V Xn2{z) J 

It is noted that the nondegenerate perturbation 
method is valid when only the two subbands with en- 
ergy Eqi and Eq2 (i.e., the lowest unperturbed subbands 



corresponding to the LH-like and SO-like hole states re- 
spectively) are accounted. Otherwise, when more sub- 
bands are included, the divergence may occur in the cal- 
culation if two energy levels Eni and En' 2 are close to 
each other, as there are terms proportional to ^ ^j^^ 
in the envelope functions "^xxniz) {x — I and s) [refer 
to Eqs. (|A12| - ET3| ]. This divergence goes to the spin- 
orbit coupling coefficients 11 and A, and finally the spin 
splitting. 

For the sake of convergence of the envelope func- 
tions and thus the spin-orbit coupling coefficients and 
spin splitting, sufficient subbands have to be considered. 
Thus, instead of the nondegenerate perturbation method, 
we apply the exact diagonalization method to obtain the 
eigen-states of Hq"'' , with sufficient basis functions con- 
structed by the two sets of functions in Eq. (|A9p [2 (20) 
of each set for Si/SiGe (Ge/SiGe) QWs]. 
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